Field Induced Spin Supersolidity in Frustrated Spin-1/2 Spin-Dimer Models 
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By means of the recently developed algorithm based on the tensor product states, the magneti- 
zation process of frustrated spin-1/2 spin-dimer models on a square lattice is investigated. Clear 
evidence of a supersolid phase over a finite regime of magnetic field is observed. Besides, critical 
fields at various field-induced transitions are determined accurately. Our work hence sheds new light 
on the search of the supersolid phase in real frustrated spin-dimer compounds. 



Over the last decade, dimerized quantum antifer- 
romagnets have attracted much attention, since rich 
field-induced quantum phases can appear in these com- 
pounds. 1 In the absence of magnetic field, the system 
exhibits spin-singlet ground state, consisting of dimers 
for closely spaced pairs of spins S = 1/2. By applying a 
large enough field such that the energy gap to spin-triplet 
excitations closes, the lowest triplet excitation starts to 
condense and the system develops a spin superfluid (SF) 
state, which supports a staggered magnetization trans- 
verse to the field direction. When the effective repulsion 
between triplets overwhelms their kinetic energy, instead 
of Bose-Einstein condensation of triplet excitations, in- 
compressible commensurate crystals of triplets with bro- 
ken translational symmetry can be stabilized, 2 which are 
signaled by magnetization plateaus. 

Recently, an even more exotic supersolid (SS) phase is 
found in the vicinity of a magnetization plateau for some 
spin-dimer models. 3 This generated an enormous interest 
in the study of quantum and thermal phase transitions 
out of these spin SS states^ The main character of 
the SS phase is that its ground states possess both solid 
and SF long range orders. It has been proposed that 
correlated hopping of triplets may play a crucial role in 
forming this novel phased At the SF-SS transition point, 
a magnetization anomaly and thereby a discontinuity in 
magnetic susceptibility appears. This indicates that the 
SF-SS transition is of second order. Moreover, approach- 
ing the edge of the magnetization plateau from the SS 
phase, the spin stiffness vanishes linearly, 4 in agreement 
with the superfluid-insulator universality class. ? 

For the spin-1/2 spin-dimer models studied previously, 
large Ising-like exchange anisotropy is necessary for the 
existence of the SS phase, 3-5 which is unrealistic for most 
magnetic systems. Under some kind of approximation, 
such an anisotropy can be considered as an effective 
interaction resulting from frustrated spin-isotropic cou- 
plings. 3,5 Thus the results for the spin-anisotropic cases 
suggest that the SS phase can appear also in frustrated 
spin-isotropic models. Nevertheless, because of underly- 
ing approximation, quantitative predictions of the needed 
frustrated coupling and magnetic field for the appearance 
of this novel phase can not be determined accurately in 
the preceding studies. To provide useful guide to the 
experimental search of the SS phase in real frustrated 



spin-dimer compounds, investigations directly within the 
frustrated spin-isotropic models are called for. 

In this paper, various field-induced quantum phase 
transitions in frustrated spin-1/2 spin-dimer Heisenberg 
models are explored by using the combined algorithm 9 
of the infinite time-evolving block decimation (iTEBD) 
method 1 ^ and the tensor renormalization group (TRG) 
approach^ 1 - In this combined algorithm, the ground 
states are assumed in the form of the tensor product state 
(TPS) or the projected entangled-pair state (PEPS). 12 
The power of the TPS/PEPS-based approach in studying 
first-order quantum phase transitions has been demon- 
strated by several groupsji 3 -"— We note that, in contrast 
to quantum Monte Carlo simulation which is plagued by 
the sign problem in studying frustrated spin systems, the 
numerical approach employed here is appropriate because 
frustration does not introduce additional difficulties to 
the TPS/PEPS-based method, 1 ^ By using the com- 
bined iTEBD and TRG method, the existence of the SS 
phase in the frustrated systems under consideration is 
firmly established, and the critical fields bounding this 
phase are determined precisely. These results provide 
further theoretical support on searching the SS phases in 
the neighborhood of magnetization plateaus observed in 
some spin dimer compounds with frustration. Moreover, 
the success in obtaining precise results for frustrated spin 
systems clearly demonstrates that the combined algo- 
rithm can be an efficient and accurate numerical tool with 
wide applications including even frustrated systems. 

We consider the following spin-1/2 bilayer Heisenberg 
Hamiltonian under a uniform magnetic field h on a square 
lattice with a strong interlayer exchange Jj_, a weaker 
intralayer coupling J, and a frustrating interlayer inter- 
action Jd, [see Fig.QJa)] 

H = J± St,l ' Si, 2 — k Si,a 

i i,a 

^ ' Sj >a • Sj a + Jd S ' Si jQ • Sj^a ■ (1) 

(%j),ot 

Here the antiferromagnetic couplings are assumed (Jj_, 
J, Jd > 0). The indices i and j denote rungs of the bi- 
layer lattice, and a = 1, 2 labels the two different layers. 
The summation runs over pairs of nearest-neighbor 
rungs. Henceforth, J± = 1 is set to be the energy unit. 
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FIG. 1: (Color online) (a) Square lattice of S = 1/2 dimers 
with an intra-dimer interaction J± and inter-dimer interac- 
tions J and Jd- (b) The construction of periodic tensor net- 
work with a 2 x 2 unit cell. 



To characterize different phases, several local order pa- 
rameters are calculated. Firstly, when systems undergo 
the dimer-SF transition, the z-component uniform mag- 
netization per site m z n = (1/2-ZV)^ (Sf 1 + Sf 2 ) begins 
to be nonzero, where N is the total number of rungs on 
the bilayer. Secondly, the SF states with nonzero in- 
plane antiferromagnetic magnetization along the x direc- 
tion can be detected by the x-component staggered mag- 
netization per site m x st = (1/2JV) £\ (Sf tl - Sf 2 ) e iC > ri . 
Using the mapping for the triplets to the semihardcore 
bosons b\ = (S^ - S+ 2 ) e lQ r ' / 'v^r tne condensate den- 
sity no of triplet excitations can be related to m* t by 
n = K&t)! 2 = 2|mg t | 2 . Finally, the checkerboard solid 
(CBS) order will be signaled by a finite value of the 
z-component staggered magnetization per site m z t = 
(1/2A0E* (Sfi + Sf t2 ) e iQl \ where Q = (tt.tt). In the 
thermodynamic limit, this quantity has a direct relation, 
\m z st \ 2 = S(Q)/2N, to the static structure factor 5(Q) = 

These quantities are evaluated under the combined 
iTEBD and TRG algorithm, 9 .^ where the ground-state 
wave function is approximated by a TPS/PEPS ansatz. 12 
Our construction of TPS for these bilayer systems on 
a square lattice is to attach a rank- five tensor d 
to each rung i and a diagonal singular value matrix 
to each bond of nearest-neighboring rungs i and 
j, as sketched in Fig. [IJb). Here s is the physical in- 
dex with s = l,---,4 for the present spin-dimer case, 
and I, r, u, d(= 1 • • • D) denote the virtual bond indices in 
four directions. The approximation can be systematically 
improved simply by increasing the bond dimension D of 
the underlying tensors. The optimized TPS/PEPS is de- 
termined through the power method via iterative projec- 
tions for a given initial state. This procedure can be con- 
sidered as a generalization of the one-dimensional iTEBD 
method™ to the two dimensional cases. To improves the 
stability the algorithm, the calculated ground state for a 
lower field h is usually employed as the initial state for 
the higher-/i case^ By means of TRG method^ the ex- 
pectation values for a TPS/PEPS ground state of very 
large systems can be calculated efficiently, where the ac- 
curacy can be systematically improved by increasing the 
TRG cutoff D cut . In the present work, we consider the 
bond dimension up to D = 5 and keep D cut > D 2 to 
ensure the accuracy of the TRG calculation. 




FIG. 2: (Color online) Values of mj, |mj t |, and no = 2jmg t | 2 
for the ground states at Jd = 0.15 (upper panel) and 0.21 
(lower panel) with J = 0.38 as functions of external field h 
for systems of size 2 7 x 2 7 with D = 4 and D cut = 16. The 
Jd independence of the saturation field h s can be observed. 



The results of the local order parameters denned above 
as functions of magnetic field h for Jd = 0.15 and 0.21 
with J = 0.38 and systems size N — 2 7 x 2 7 are shown 
in Fig. [5J Here we take the bond dimension D = 4 and 
the TRG cutoff D cut = 16. We note that results for 
D = 3, 5 are very similar to those for D — 4 (see below). 
Upon increasing the field, the ground state of the model 
in Eq. ([I]) goes through a succession of distinct quantum 
phases. For the case of weak frustration (Jd = 0.15), the 
ground state evolves first from the dimer phase (m^, no = 
0) to the SF phase (m*, no ^ 0), and finally to the fully- 
polarized (FP) phase (to* — 1/2 and n = 0). Such a 
behavior is consistent with the results of the unfrustrated 
(J d = 0) case3 

For the larger- Jd case (Jd = 0.21), apart from afore- 
mentioned phases, a magnetization plateau and a pecu- 
liar SS phase with both the CBS and the SF orders are 
found. The general feature of our results for Jd = 0.21 
is quite similar to what have been discovered in the pre- 
vious works for the spin anisotropic systems^— Within 
the regime of magnetization plateau, we find that the uni- 
form magnetization is m„ ~ 0.25 together with a nonzero 
CBS order parameter |mg t | ~ 0.23. Thus this plateau 
state describes the CBS phase with long-range diagonal 
order at wave vector Q = (tt, it). The observed CBS 
phase can be understood as formed by half-filling the sys- 
tem with triplet excitations so that rnj; = 0.25. Its stabil- 
ity comes from strong effective repulsion between near- 
est neighbor triplets generated by interdimer frustrated 
couplings J and Jd- We note that, as in the anisotropic 
cases studied previously^— the calculated CBS order pa- 
rameter \m z t \ is also slightly smaller than its classical 
value |TOg t | = 1/4, which is caused by quantum fluctua- 
tions therein. At higher fields, the melting of the CBS is 
found to be of first order and the SF phase reappears. At 
this transition, the magnetization changes abruptly from 
~ 0.25 to ~ 0.35, and the condensate fraction no 
jumps from zero to no ~ 0.15. On the other hand, when 
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FIG. 3: (Color online) Dependence of the calculated order 
parameters on D cut (upper panels) and D (lower panels): (a) 
h = 1.2 and (b) 1.3 with D = 4; (c) h = 1.2 and (d) 1.3 with 
-Dcut = 16 for D < 4 and D cut = 25 for D = 5. 
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FIG. 4: (Color online) Details of local order parameters 
around (a) the SS phase, (b) the dimer-SF and (c) the SF-FP 
transition points with enlarged scale. The parameters are the 
same as those for the lower panel of Fig. [2] The inset in panel 
(a) shows the critical fields h C 2 (x) and /i C 3 (*) for various D 
(Dcut = 16 for D < 4 and £> cut = 25 for D = 5). 



decreasing the field from the plateau, a SS phase emerges 
in the present frustrated spin-dimer system, where both 
the CBS order parameter m| t and the Bose condensate 
density no are nonzero. At even lower fields, the CBS or- 
der disappears and a standard SF phase with only no 
is recovered. 

As a check on our prescription of TRG cutoff D cut 
and bond dimension D, the dependence of the evaluated 
order parameters on them is shown in Fig. [3] Here we 
consider two typical values of h = 1.2 and 1.3, which are 
corresponding to the SF and the SS phases respectively. 
It is expected that the ground states in these two cases 
are more entangled than those in the dimer or the CBS 
phases. Consequently reliable results can be reached only 
when D and D cut are large enough. As seen from Fig. [31 
there is almost no dependence on D cut from D cut = 10 
to 16 for D = 4, and the dependence on D is small. (The 
variance of each order parameter is less than 0.025 even 
in the SS phase with h = 1.3.) These results indicate 
that the bond dimension and TRG cutoff used in this 
work is large enough to provide accurate findings for the 
present case. 

To explore the behaviors around the transition points 
to the SS phase, results around this phase with enlarged 
scale are shown in Fig. Ufa) . We find that the transition 
between the SF and the SS phases occurs at the critical 
field h C 2 ~ 1.26, at which the CBS order mj. begins to 
vanish. As observed in Fig. Eta) (also in Fig. this 
transition is also signaled by a kink in the magnetization 
curve, which implies a discontinuity in magnetic suscepti- 
bility x = dm^/dh. The field-induced transition between 
the SS and the CBS phases occurs at another critical field 
h C 3 1.34, at which the condensate fraction no begins 
to vanish. Moreover, we observe that no vanishes lin- 
early as h approaches to h c ^ from below. This behavior 
is in agreement with the SF-insulator universality class,® 
since the density of the bosonic excitations (triplet holes) 
should be small around h c ^ and the condensate fraction 



roughly equals to the spin stiffness in this dilute-boson 
limit. In the inset of Fig.QJa), the dependence of the crit- 
ical fields h C 2 and h c % on bond dimension D is presented. 
The values show only a minor dependence on D, which 
again supports the validity of the present approach. 

Now we turn our attention to the cases near the dimcr- 
SF transition at the critical field h c \ and around the sat- 
uration field h s . In both cases, some analytical results 
are available. As same as the unfrustrated (Jd = 0) 
caseji^ the instability of the dimer phase is triggered 
by energy gap closing of the bosonic one-triplet exci- 
tation above the dimer state. On the other hand, the 
FP state becomes unstable when the lowest bosonic ex- 
cited state with a single spin flip is degenerate with the 
FP state. Thus both transitions are expected to be of 
second order and in the universality class of the dilute 
Bose gas quantum-critical points By using the expres- 
sion of the excitation spectrum of the one-triplet states 
up to the third-order perturbation expansion, 2 — the criti- 
cal field h c i of the dimer-SF transition is achieved, h c \ ~ 
1 — 2(J — Jd) — | J( J— Jd) 2 ■ To determine the saturation 
field h s , one notes that the states with a single spin flip 
are the eigenstates of our model in Eq. (JXJ) and their en- 
ergy eigenvalues can be calculated exactly. Thus the ex- 
act expression for h s can be derived. We find that, when 
Jd < J and Jd < 1/4, the saturation field h s — 1 + 4 J, 
which is independent of the magnitude of the frustrated 
coupling Jd- A closer look around the transitions at h c ± 
and h s is shown in Figs. EJb) and (c). In agreement 
with the SF-insulator universality class^ we find that 
the condensate fraction no vanishes linearly also around 
these two field-induced transition points. The dimer-SF 
transition occurs at h c i ~ 0.61, where the uniform mag- 
netization m u and the condensate fraction no start to be 
nonzero. Above the saturation field h s w 2.52, all spins 
are polarized in the z direction and the uniform magne- 
tization saturates at the value of ?n u = 0.5. Our findings 
of h c i and h s agree well with the values (h c \ ~ 0.64 and 
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FIG. 5: (Color online) Local order parameters at Jd = 0.23 
(upper panel) and 0.27 (lower panel) as functions of external 
field h. Other parameters are the same as those for Fig. [21 In- 
sets of both panels show the details of local order parameters 
around their SS phases with enlarged scales. 

h s = 2.52) given by the analytic formulas with J = 0.38 
and Jd = 0.21. Moreover, as seen from Fig. [21 the cal- 
culated values of the saturation field h s for the cases of 
Jd = 0.15 and 0.21 are identical. Thus the independence 
of the saturation field on the frustrated coupling Jd is 
indeed preserved in our approach. These observations 
further substantiate that the employed method can pro- 
vide reliable results even for the frustrated spin systems. 

As shown in Fig. [21 a novel SS phase can exhibit as 
ground states of frustrated spin-dimer models when the 
frustrated coupling Jd becomes large enough. One may 
wonder what the fate of the SS phase is if Jd increases 



further. In Fig. [5l two values of even larger Jd {Jd = 
0.23 and 0.27) are considered. We find that the width 
in h of the SS state slightly decreases from h c s — h C 2 ~ 
0.08 for J d = 0.21 (see Fig. ij down to about 0.06 for 
Jd = 0.27 (See Fig. [5]). Nevertheless, as seen from Fig. [SJ 
the regime for the plateau state expands upon increasing 
Jd- Moreover, the width in h of the SF phase shrinks 
and the magnitude of no in this state decreases as Jd 
increases. These outcomes indicate that the SS state may 
eventually disappear for large enough frustration. That 
is, the SS phase may be stabilized only within a limited 
region of the h- Jd phase diagram. Therefore, carefully 
tuning system parameters into the suggested parameter 
regime are necessary to uncover experimentally this novel 
phase in real frustrated spin-dimer compounds. 

To summarize, the field-induced quantum phase tran- 
sitions in frustrated spin-1/2 spin-dimer models are stud- 
ied under the combined iTEBD and TRG algorithm!*^ 
Clear evidence of a SS phase over a finite regime of mag- 
netic field is provided. Critical fields at various field- 
induced transitions are evaluated accurately by using 
merely moderate bond dimension D. Besides providing 
quantitative predictions of the needed frustrated coupling 
and magnetic field for the appearance of the SS phase, 
the present work also demonstrates clearly the potential 
in applying the current formalism even to frustrated spin 
systems. 
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